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Abstract. The discovery of the fractional quantum Hall effect in GaAs-based semiconductor 
devices has lead to new advances in condensed matter physics, in particular the possibility for 
exotic, topological phases of matter that possess fractional, and even non-Abelian, statistics 
of quasiparticles. One of the main limitations of the experimental systems based on GaAs 
has been the lack of tunability of the effective interactions between two-dimensional electrons, 
which made it difficult to stabilize some of the more fragile states, or induce phase transitions 
in a controlled manner. Here we review the recent studies that have explored the effects of 
tunability of the interactions offered by alternative two-dimensional systems, characterized by 
non-trivial Berry phases and including graphene, bilayer graphene and topological insulators. 
The tunability in these systems is achieved via external fields that change the mass gap, or by 
screening via dielectric plate in the vicinity of the device. Our study points to a number of 
different ways to manipulate the effective interactions, and engineer phase transitions between 
quantum Hall liquids and compressible states in a controlled manner. 



1. Introduction 

Strongly correlated phases of electrons confined to move in the plane, subject to a perpendicular 
magnetic field, have attracted significant attention since the discovery of fractionally quantized 
Hall conductivity [1]. The profound role of topology in this extreme quantum limit leads to the 
presence of quasiparticles that carry a fraction of electron charge [2], and fractional (Abelian, 
or possibly non-Abelian) statistics [21 O HI El \7\. The prospect of excitations possessing 
non-Abelian statistics has motivated different schemes for topological quantum computation [8] 
based on these systems. 

These remarkable phenomena occur in the fractional quantum Hall (FQH) regime, when the 
number of electrons, N e , is comparable to the number of magnetic flux quanta iV$ through 
the two-dimensional electron system (2DES). Correlated FQH liquid states appear at certain 
partial filling v = N e /N$ of the "active" Landau level (LL). In traditional semiconductor 
heterostructures, the physics of a partially-filled n = 1 LL differs significantly from that of 
n = LL, due to the node in the single-particle wavefunction [9j. As a consequence, the 
hierarchy /composite fermion states [101 111], ubiquitous in the lowest Landau level (LLL), are 
significantly weakened in n = 1 LL, and some of the more exotic states, such as the Read-Rezayi 
(RR) parafermion states [12], are likely to be favored. A number of studies have focused on the 
simplest, k = 2 non-Abelian member of the RR sequence - the Moore-Read (MR) "Pfaffian" 
state [5], believed to describe the FQH plateau at v = 5/2 [13]. Quasiparticles of the MR and 



higher RR states obey the non-Abelian statistics [5] which is of interest for topological quantum 
computation [8]. 

One of the main disadvantages of the GaAs-based devices is that their 2DES is buried inside 
a larger, three-dimensional structure. This unfortunate fact fixes the effective interactions at 
values that are often not optimal for some of the most interesting FQH states, including the RR 
series. For example, the MR state is found to lie very close to the boundary with a compressible 
phase |14[ [15] . Another problem stems from strong dielectric screening and finite well- width [TB] 
in GaAs, which weaken the electron-electron interactions and make FQH states fragile. This has 
been a major obstacle in the studies of the non-Abelian states, which could only be observed in 
ultra-high-mobility samples [12]. Thus, it is desirable to find an alternative high-mobility 2DES 
with strong effective Coulomb interactions that are adjustable in a broad range. 

Recently, a new class of such high-mobility 2DES which host chiral excitations with non- 
trivial Berry phases, has been discovered. These "chiral" materials include graphene and 
bilayer graphene [T7], and, more recently, topological insulators |18j . as well as certain quantum 
wells [19] . The chiral nature of the quasiparticles gives rise to new electronic properties, including 
an unusual LL sequence, the anomalous Hall effect, and suppression of weak localization |17j . 
When these chiral materials are subject to a perpendicular magnetic field, the kinetic energy 
quenches into discrete Landau levels, similar to the usual semiconductors with non-chiral carriers. 
However, due to the chiral band structure [20J and the fact that the surface of these materials is 
exposed [21], they offer new possibilities to tune the effective interactions and explore strongly 
correlated phases. 

In this paper we review two practical ways |21[ [20] of tuning the effective interactions in chiral 
2DES, and the effect this has on the FQH states. One of the attractive features of the chiral 
2DESs is that they allow for a more robust realization of certain FQH states in several LLs, 
as opposed to a single LL in GaAs. Additional insights can be obtained by driving transitions 
between the incompressible (FQH) states and the compressible, stripe and bubble, phases. Such 
transitions can be implemented in chiral and massive 2DESs by varying the external field. 
Overall, the tunability of the chiral 2DES allows one to explore a larger region of the effective 
interactions than has been achieved in GaAs. 



2. Model 

We begin this Section by providing a brief overview of the one-body problem in a magnetic field, 
as well as the summary of the standard recipe of exact diagonalization applied to quantum Hall 
systems. 

Consider a single electron moving in a plane, subject to a perpendicular magnetic field 
Bz = V x A(r). The corresponding Hamiltonian can be written in a "covariant" notation 
as 

K = ^-g ab n a n b , (1) 
2m 

where ir a = p a — -A a (r) (a = x,y) represents the dynamical momentum, and g is the band 
mass tensor. In the usual isotropic case {g = I), we can obtain the single-particle energies 
(Landau levels) by choosing, for example, the symmetric gauge A x = By/2,A y = —Bx/2. In 
this case, the dynamical momenta become tt x = —ih-^ + ^fry and -K y = —ih-^j — -^rx, in terms 

of the magnetic length £b = yh/eB. Hamiltonian can be transformed into a diagonal form 
K = ^ (a) a + \) with the help of "ladder" operators 

a oc ir x + iiTy, (2) 
oc tt x — i-Ky. (3) 



However, for each value of era, there remains a residual degeneracy equal to the number 
of magnetic flux quanta, iV$. This degeneracy is resolved by another pair of operators 
b, that commute with a, and depend on the guiding center coordinates of the electron, 
R a = r a - 

^fi'^b^B- Operators b^ create the (unnormalized) single particle eigenstates of the 
lowest, n = Landau level (LLL), 

Mz)Kz l e- z * z ^B, (4) 

with z = x + iy being the complex coordinate of an electron in the plane (and z* denoting its 
complex-conjugate). The quantum number I is an eigenvalue of the angular momentum, and 
the single particle states are localized along concentric rings around the origin. Single-particle 
states in higher LLs are constructed by the repeated action of cr, thus completing the solution 
of the one-body problem in a rotationally-invariant case. 

2.1. Many-body problem: Exact diagonalization 

An interacting problem of a finite number of electrons can be studied numerically, using exact 
diagonalization which has been remarkably successful in unravelling many of the essential 
physical properties of FQH systems, for systems as small as 10 electrons [21 [TO]. A crucial 
approximation that makes exact diagonalization practical is to neglect the excitations between 
LLs and restrict the Hilbert space to a single, "active" LL, which has a degeneracy of N$. 
This approximation is physically justified in high magnetic fields where such excitations indeed 
become very costly. One can then construct a many-body Hilbert space, consisting of Slater 
determinants I2, ■ ■ ■) built from the states (|3J), using a suitable choice of boundary condition. 
In the study of the FQH effect, two kinds of surfaces are available that preserve the translational 
invariance of an infinite 2DES: sphere |10j and torus [22]. The two choices of boundary conditions 
illustrate the specific features of a many-body FQH state under investigation: on a sphere, the 
FQH state couples to the curvature of the manifold, which is characterized by the topological 
number called shift [23] . The shift produces a small offset between N$ = N/v and the magnetic 
monopole, whose strength is denoted by 2S, placed in the center of a sphere. As a consequence, 
different FQH states "live" in different Hilbert spaces and in principle can be directly compared 
only after extrapolation to the thermodynamic limit. On the other hand, the flat surface of a 
torus leads to a unique definition of iV$ for given N and v, such that different candidate states 
describing the same filling v are all realized in the same Hilbert space. The caveats of this 
geometry are the additional geometric parameters, the angle and aspect ratio of the torus; the 
Hamiltonian of a finite-size system depends on these parameters and their specific values favor 
one FQH phase over the others. Thus, an analysis is slightly more involved but the gain is that 
the ground-state degeneracy can be used to identify topologically-ordered states. In contrast, 
on the sphere FQH states always appear as non-degenerate, zero angular momentum ground 
states. 

Once the Hilbert space is defined and properly adapted to the symmetry group/boundary 
condition, the remaining task is to represent the Hamiltonian and diagonalize the corresponding 
matrix using a variant of Lanczos algorithm [24J (or Jacobi/Householder routine [25] for very 
small systems). The advantage of symmetry, in particular rotational invariance, is that Wigner- 
Eckart theorem applies, and leads to the definition of the so-called Haldane pseudopotentials [9|. 
Any two-body interaction, such as Coulomb potential projected to a single LL, can be written in 
the usual second-quantized form H = X^{m 4 } ^mim2m3m 4 Cm 1 Cm 2 c m4 c m3 . Furthermore, assuming 
rotational invariance, the matrix element factorizes into 

2S L 

v .-Vv, V r SSL r SSL (*\\ 

L=0 M=-L 



i.e. into a product of a geometrical factor involving the Clebsch-Gordan coefficients C^^ m , 
and the Haldane pseudopotential Vl • The finite set of Haldane pseudopotentials contains all the 
information about the interaction projected to a single LL and thus determines all the physical 
properties of a FQH system. Haldane pseudopotentials can be conveniently evaluated from the 
Fourier transform of the interaction, V(q). 

2.2. Modification of the interaction due to the spinor nature of the wavefunction in massive 
Dirac materials 

Once the Hilbert space is defined and symmetry-adapted, the Coulomb Hamiltonian can be 
represented in the matrix form. In the LLL, the interaction is simply given by V(q) = 1/q, which 
is the Fourier transform of the Coulomb potential in two dimensions. However, if we consider the 
case when the chemical potential is such that LLs from 1 through n — 1 are completely filled and 
"inert" , while the nth LL is partially filled and "active" , the effective interaction becomes slightly 
modified and assumes the form 1/q x \F n {q)\ 2 , where the form factor [9] F n (q) = L n (q 2 £ B /2) 
can be evaluated using the algebra of operators a, a), Eq. ([2]) - see, e.g., Ref. [TTJ (L n is the 
nth Laguerre polynomial). Similarly, if the 2DES has a finite extent in the z-direction, the 
interaction becomes modified via an appropriate form-factor that results from integrating out 
the z-component of the wavefunction. 

The modification of the effective interaction may also result from the spinor nature of the 
wavefunction |26| . Such wavefunctions arise naturally in a family of 2D materials that have non- 
trivial Berry phases. One such material is monolayer graphene, a high-mobility atomically thick 
2DES [T7j, where recently several FQH states of the type v = m/3 have been discovered |27j . 
A closely related material, bilayer graphene |17| has similarly high mobility, and exhibits 
interaction-induced quantum Hall states at integer filling factors at low magnetic fields [28] . 
Graphene and its bilayer are characterized by Berry phase tt (graphene-like) and 2n (bilayer- 
graphene-like with an energy gap), respectively. More explicitly, we refer to the Dirac materials 
as those described by a family of 2 x 2 Hamiltonians with Berry phase tt and 2tt |20] . The case of 
7r-carriers is realized in graphene, topological insulators, special quantum wells [TS]; the case of 
27r-carriers occurs in bilayer graphene, where the energy gap can be controlled in a wide range 
by a perpendicular electric field [17 . Details on how the effective 2x2 Hamiltonian can be 
derived from a more complete tight-binding model can be found in the review |29j . 

For 7r carriers, the single-particle wavefunctions are given by spinors [20] ip n = 
(cos n 0| n |_i, sin n 0| n |), where <f) n is the wave function of the nth non-relativistic LL (n = 
±1, ±2, . . .), and parameter 9 depends on the ratio A/(Hvo/£b), where A is the mass gap, and 
vq the Fermi velocity. We use the notation (<77r, n) to denote the nth LL for air carriers (a = 1,2). 
As a consequence of the spinor wavefunction, the effective form factor [9] F£(q) that describes 
the interaction projected to a (tt, n) LL is given by 

F-(q) = cos 2 9L ]nhl (t^>j + sin 2 9L ]n] , (6) 

where for simplicity we omitted the index of 6. The form-factor is a mixture of the (|n| — l)th 
and |n|th LL form-factors in a non-relativistic 2DES with parabolic dispersion. At 9 = vr/4, the 
above equation reduces to the form-factor of graphene [26], however, by varying /S./{Kvq/£b) in 
the experimentally accessible range [20], one can realize any value of £ (0;7r/2). 

Similarly, for carriers with Berry phase 27r, the single-particle wavefunctions are ip n — 
(cos n (f>\ n ii, sin O n (f>\ n \ + i), and the form- factor is a mixture of standard (|n| — l)th and (|n| + l)th 
form- factors, 

F^(q) = cos 2 9L lnhl (^) + sin 2 0L |n|+1 (^) . (7) 




Figure 1. (Color online) Two ways of modifying the effective interaction V{q) in the FQH 
regime: (a) via dielectric screening, for several values of a given in the inset (using Eq.([8])) for 
d = £b, (b) via chiral band structure, for several values of cos 2 6 given in the inset (in (tt, 1) LL). 
Because the interaction V(q) is isotropic, it is plotted as a function of q x only. The comparison 
of the corresponding Haldane pseudopotentials for cases (a) and (b) can be found in Ref. [30J. 



The tunable form of the effective interactions Eqs. (|6l7p provides a way to engineer transitions 
between strongly correlated phases in situ by changing the field, as illustrated in Fig. Hfb). 



2.3. Modification of the effective interaction due to dielectric screening 

Up to this point, we have discussed the change in the effective interaction V(q)\F(q)\ 2 that 
resulted from a modification of the single-particle wavefunctions, and therefore F(q). The 
method, proposed in Ref. [2T], allows one to directly change V(q). Chiral materials, such 
as graphene, are often exposed to the environment, which allows for dielectric material to be 
deposited on top of them. We consider a setup where graphene sample is situated in a dielectric 
medium with permittivity ei, and a semi-infinite dielectric plate with permittivity e<i ^ e\ is 
placed at a distance d/2 away from the graphene sheet. The effective interactions between 
electrons in graphene change due to the surface charges induced at the boundary between 
dielectrics [21]: 

e 2 e 2 ei — €2 

V{r) = h q , where a = . (8) 

t\r eiVr 2 + d 2 ei + £2 

The ratio d/is controls the effective interactions within a partially filled LL. However, the 
overall energy scale is also modified and this has an impact on the magnitude of the excitation 
gap. The gap should be multiplied by a factor ecaAsAi if comparison is to be made with GaAs 
2DES. Again, an important advantage of this setup is that the interactions can be tuned in 
situ by varying the magnetic field B, which modifies the ratio d/ls- The consequences for the 
many-body system are illustrated in Fig. [Ha). 



3. Results 

We now present some exact-diagonalization results for the interaction models introduced in Seci2J 
The focus is on the physical features resulting from the variation of the effective interaction; 
as stated in the introduction, we follow the traditional assumption of neglecting the excitations 
between different LLs, which is expected to be quite reasonable in GaAs, but perhaps less so in 
chiral materials. Much of previous theoretical work on chiral materials has been restricted to 
graphene, exploring the consequences of the four-fold LL degeneracy (valley and spin) that leads 
to new SU(2) and SU(4)-symmetric fractional states [31J. Instead, we focus on the high-field 
limit, neglecting the multicomponent degrees of freedom, and examine the effects originating 
from the interplay of the Coulomb interaction and band structure. We consider a family of 



Figure 2. (Color online) Phase transition between the Laughlin and bubble phase in (2tt, 1) LL. 
(a) Low-lying energy spectrum on the torus as a function of 9. (b) Projected structure factor 
Sl as a function of angular momentum L on the sphere and parameter 9. Large ground-state 
degeneracy for small cos 2 9, as well as the sharp peak in Sl, correspond to the bubble phase. 

band structures introduced in Ref. [20], which describe a number of high-mobility materials, 
including graphene with massive carriers (mass is generated either spontaneously, or as a result 
of sublattice symmetry breaking [17]). topological insulators [18], bilayer- graphene [17] . trilayer 
graphene |32j, and similar materials. Pristine graphene, which hosts massless Dirac-like fermions, 
is contained in this model as a particular case. Note that we explore fairly large variations of the 
effective interactions among fully polarized electrons; more subtle effects due to LL mixing [33 
and multicomponent degrees of freedom are left for future study. 

We show the results for both spherical and torus geometry. A number of useful insights 
can be inferred from the study of the energy spectrum; in addition, we use overlap calculations 
to compare an exact, many-body ground state ^exact) with a numerical representation of a 
trial wavefunction, Atrial- The overlap is defined as a scalar product between two normalized 
vectors, O = |( Atrial Inexact) | • If O is consistently close to unity for a number of system sizes 
considered, we consider the trial wavefunction to be a faithful description of a FQH phase. From 
the knowledge of a ground-state wavefunction, we also evaluate the (projected) static structure 
factor [35]. Sharp peaks in the structure factor indicate the onset of compressible phases [Ml [36]. 

The simplest example of the phase transition can be realized between the Laughlin state and 
a bubble phase in (2ir, 1) LL of the chiral materials. In Fig. [2]^a) we show the low- lying part 
of the energy spectrum obtained by exact diagonalization in the torus geometry, as a function 
of parameter 9 introduced in Sec. [2j For values cos 2 9 > 0.2, the ground state of the system is 
non-degenerate and has a large overlap with the Laughlin wavefunction. When cos 2 9 is less than 
0.2, the effective interaction resembles that of n = 2 non-relativistic LL, which favors the bubble 
phase at filling v = 1/3. This is manifested by a large ground-state degeneracy developing in 
the spectrum [36J. Similarly, the evolution of the projected structure factor Sl on the sphere as 
a function of 9 develops a discontinuity around cos 2 9 ~ 0.2 and a sharp peak below that value, 
signaling the compressibility of the ground state. 

Previously [20], it was emphasized that the tunability of the interactions in chiral materials 
can lead to the more stable non-Abelian states. One such example of considerable experimental 
interest is the Moore-Read Pfaffian state [5]. In Fig. [3] we show the energy spectrum for A = 12 
electrons on a torus in (tt, 1) LL. At cos 2 9 = 0, the interaction reduces to n = 1 non-relativistic 
Coulomb interaction, which leads to a fragile Pfaffian state. As shown in the inset of Fig. [3l 
the overlap of the exact ground state with the Moore-Read wavefunction slightly increases as 
cos 2 9 is increased from zero. However, due to the particle-hole symmetry which is present on 
the torus at filling v = 1/2, a better indicator is to compare the exact ground state with a 
particle-hole symmetrized version of the Moore- Read wavefunction [15]. Indeed, in this case 
the overlap rises to nearly unity right before the transition at cos 2 9 ~ 0.14. At this point, the 
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Figure 3. (Color online) Enhanced Pfaffian state and the transition to the CFL. Energy 
spectrum of N = 12 electrons on the torus in (tt, 1) LL is plotted as a function of 6. Inset shows 
the ground-state overlap with the Moore-Read wavefunction and its particle-hole symmetrized 
version. 

level crossing occurs, and a composite Fermi liquid (CFL) state becomes the ground state of the 
system. We emphasize that although the level crossing would suggest a first-order transition, 
we cannot reliably rule out a possibility for the second-order transition in the thermodynamic 
limit. Note that, in addition to the considerable increase in the overlap with the Moore-Read 
state, also the excitation gap increases towards the transition to the CFL state, adding to the 
overall stability of the Pfaffian. As shown in Ref. [20J, Pfaffian correlations can also be realized 
in (tt, 2) LL, and other states such as k = 3 RR state can be realized in several LLs [30j . 

Phase transitions can also be studied without reference to a particular wavefunction by 
evaluating the excitation gap of the system, Fig. 0] To be more general, we consider a model 
that involves a superposition of n = 0, n = 1 and n = 2 LL non-relativistic form factors [20j: 

F n (q) = cos 2 0L w (^) +sin 2 #cos 2 </>L w+1 + sin 2 9 sin 2 (ftL ]n]+2 (^) . (9) 

This model generalizes the case of (tt, 1), (vr,2) and (2tt, 1) LLs studied earlier, and naturally 
arises in systems such as trilayer graphene [32]. In Fig. 2] we plot the charge gap at v = 1/3 and 
v = 1/2 as a function of and (ft. Because the calculation is performed in the spherical geometry, 
we fix the shift to be —3 (corresponding to the Laughin state at v = 1/3 and the Moore-Read 
state at v = 1/2). Along certain lines, indicated by arrows, the phase diagram reduces to one 
of the cases mentioned above. Note that proper finite-size scaling needs to be performed in 
order to get the correct values for the gap; however, it was previously found [21, 20J that this 
rigorous analysis produces values that are roughly in agreement with the ones shown in Fig. [5J 
At v = 1/3 the Laughlin state is realized for a wide range of 9 and (ft. The maximum of the 
gap occurs for the effective interaction that is a mixture of n = and n = 1 non-relativistic 
form-factors. For cos 2 6 < 0.2 the excitation gap is significantly reduced and a bubble phase 
is realized. The point of transition to the bubble phase is insensitive to the value of (ft, unless 
cos 2 (ft ~ 1. On the other hand, at v = 1/2 the behavior of the gap suggests the existence of 
three phases. On the right side (large cos 2 #), the gaps are close to zero - this is the region 
corresponding to the CFL phase. In this parameter regime, the ground state of the system 
also has a different shift from the one chosen in Fig. HI On the left side, one can discern two 
distinct branches - the middle one corresponds to the Moore-Read state because the overlap 
with the Pfaffian wavefunction is high in this region [20] . Note that this region is fairly large and 
encompasses the effective interactions which are quite different from the n = 1 non-relativistic 
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Figure 4. (Color online) Charge gap (in units of e 2 /eis) at z/ = 1/3 and j/ = 1/2 as a function of 
the band-structure parameters 9 and 4>. The calculation is performed in the spherical geometry, 
for TV = 10 particles at v = 1/3 and N = 14 particles at v = 1/2, with a fixed shift of —3. 

Coulomb interaction, known to give rise to the Pfaffian state. On the far left there is a stripe 
phase; the excitation gap is seemingly large in this region, however this is only an artefact of the 
spherical geometry that cannot adequately accommodate the broken-symmetry state such as 
the stripe. We expect this gap to be significantly reduced as one approaches the thermodynamic 
limit. Note that the gaps of the incompressible states can be further increased using a dielectric 
setup proposed in Ref. [21] . 

4. Conclusion 

In this paper we have discussed two practical ways of tuning the effective interactions in fractional 
quantum Hall systems. We emphasized an important role of chiral 2DESs where the interactions 
can be varied in a broader range than in GaAs, and more robust non-Abelian states are likely 
to be found. Moreover, these non-Abelian states are expected to exist in the new regimes of 
the effective interactions that are significantly different from the ones in GaAs. Finally, their 
excitation gaps can be tuned, and quantum phase transitions can be engineered in situ by 
varying the magnetic field. 
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